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o. The absorption spectrum of CO2 in the wavelength range 120 nm — 160 nm is an- 



< 



alyzed by means of quantum mechanical calculations performed using vibronically 

^^ coupled PESs of five singlet valence electronic states and the coordinate dependent 

O, transition dipole moment vectors. The thermally averaged spectrum, calculated for 

G T = 190 K via Boltzmann averaging of optical transitions from many initial rota- 

(~| tional states, accurtely reproduces the experimental spectral envelope, consisting of a 

q 

cy^ low and a high energy band, the positions of the absorption maxima, their FWHMs, 

'^ peak intensities, and frequencies of diffuse structures in each band. Contributions 

,J^ of the vibronic interactions due to Renner- Teller coupling, conical intersections, and 

'""' the Herzberg- Teller effect are isolated and the calculated bands are assigned in terms 

<^ of adiabatic electronic states. Finally, diffuse structures in the calculated bands are 

\^ vibronically assigned using wave functions of the underlying resonance states. It 

CN is demonstrated that the main progressions in the high energy band correspond to 

(^ consecutive excitations of the pseudorotational motion along the closed loop of the 

m 

^—1 CI seam, and progressions differ in the number of nodes along the radial mode per- 

.J^ pendicular to the closed seam. Irregularity of the diffuse peaks in the low energy 

5^ band is interpreted as a manifestation of the carbene-type 'cyclic' OCO minimum. 

I. INTRODUCTION 

This paper describes the results of an ab initio quantum mechanical study of the ab- 
sorption spectrum of carbon dioxide photoexcited with the ultraviolet (UV) light between 
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120 nin and 160 nin and offers vibronic assignments of the diffuse absorption bands. A brief 
account of this work has already been published-l^ The potential energy surfaces (PESs) used 
in these calculations are described in the preceding paper (termed 'paper r)P 

Photodissociation of carbon dioxide is of considerable importance for atmospheric and 
planetary chemistry. CO2 is the second common trace gas in the Earth atmosphere and 
one of the main products of the fossil fuel burning. Its emissions continue to grow^l despite 
considerable effort to mitigate climate change.^EJ UV photodissociation destroys CO2 with 
a unit quantum yield, and in the upper atmosphere, this destruction channel is dominant. 
Photolysis, shifted into the UV/visible range via preheating of CO2, has been used in solar 
based schemes of conversion of atmospheric CO2 into fuel.'^ UV photochemistry of carbon 
dioxide controls chemical processes in the upper layers of the CO2 based atmospheres of 
Mars and Venus.'^ Accurate predictions of the low temperature UV absorption cross sections 
are required to improve the existing photochemical models of their atmospheres, and are 
also relevant for the models of the atmosphere of Titan,!^ which contains CO2 as a minor 
constituent. Finally, photodissociation shares its activated complex — the highly electron- 
ically and vibrationally excited CO2 — with another reaction, in which C + O2 collisions 
produce carbon monoxide, the second most abundant molecule in interstellar clouds.'^ 

Absorption spectrum of CO2 in the UV range is studied experimentally since the begin- 
ning of the 20th century.'^^"^ Figure CTa) gives an overview between 220 nm and 105 nm. 
Carbon dioxide is transparent down to about 220 nm (6.20 eV). Below 220 nm, the absorp- 
tion starts with a very weak progression of narrow lines .'^^'^^'^ At 112.1 nm (11.08 eV), a 
strong absorption band is observed, which is assigned^^*^ to the optically allowed transition 
to the Rydberg state ^S+; this band is the first in a series of Rydberg transtions at shorter 
wavelengths [not shown in Fig. flla)]. Between 160nm and 120nm (7.75eV — 10.33eV), 
two relatively weak absorption bands are located, with the stronger one centered at 133 nm 
(9.30 eV) and the weaker one at 148 nm (8.40 eV). Magnified view of this spectral interval is 
given in Fig. W(h)- Each band exhibits irregular diffuse structures superimposed on broad 
structureless background. The vibronic origin of the diffuse structures and the dissociation 
mechanisms in these two bands are subjects of this paper. 

Previous analysis established key electronic states in the observed spectrumJ ^^ * ^^ * ^^" ^ 
These are the singlet valence states 2, 3^A' and 1, 2, 3^A", properties and topographies of the 
global PESs of which are discussed in paper I. It is customary to give electronic assignment 



using the D^oh labels, appropriate for linear CO2, namely l^S^ , l^H^, and l^A^, with the low 
energy band at 148 nm assigned to the 1^A„ state, and the high energy 133 nm band to the 
l^IIg state.'^^'^^'^ However, the interstate vibronic couplings calculated in paper I are strong 
rendering such 'zeroth order' diabatic labels somewhat problematic. Vibrational progressions 
are isolated using solely line spacings. For the stronger lines near 133 nm, one bending or 
two stretching progressions are offered as alternative tentative assignments.'^S'^Il xhe weaker 
band near 148 nm is presumed to involve strong bending excitations. This assignment is 
based on the analysis of DixorP^ who characterized the upper state in the 148 nm transition 
as bent (two-dimensional contour maps of the lowest two bent excited states 2M' and iM" 
are shown in Figs. 9 and 13 of paper I). However, the intervals between diffuse peaks are 
strongly perturbed and rather irregular in both absorption bands: A unequivocal assignment 
is difficult without a convincing spectroscopic model or accurate dynamics calculations. 

Irregularities in the peak spacings and intensity distributions are due to strong vibronic 
interactions, rooted in the symmetry properties of the excited valence states. In the near 
linear Franck-Condon (FC) region, vibronic effects can be grouped into three classes: 

(1) The rotoelectronic Renner- Teller (RT) effect involves A' and A" components of the 
orbitally degenerate l^H^ or 1^A„ state. Figure p[a) illustrates how the degeneracy is lifted 
as CO2 bends. In the shown cut, the upper two states stem from l^H^, the lower two 
stem from 1^A„, and the transitions l^Hg ^ -^^^t ^md l^A^ ^ -^^^t "^^^ ^le nominally 
classified as linear-linear and linear-bent, respectively. In either case, the A' /A" interaction 
~ AQ/ sin^ Ooco is proportional to the projections A and Q of the electronic (L) and the 
total angular momentum (J = A^co2 + L) on the molecular axis and, for Q ^ 0, diverges as 
«oco — > X80°.112128H30] jf J jg normal to the CO2 axis and fi = 0, the RT coupling is quenched. 
Interaction between A' and A" states in this case is mediated only by the Coriolis coupling 
which allows transitions from f2 = to fi = 1 block in which RT coupling is operative. 

(2) The symmetry allowed conical intersections (CIs) occur between states 1^A„ and l^H^; 
two avoided crossings in Fig. 0(b) near 2.2 ao and 2.9 oq illustrate positions of the CIs. 
Calculations detect CIs for two A' states, 2M' and 3M', and two A" states, 1, 2^A" or 2, 3M". 
In both symmetry blocks, the degeneracies are lifted linearly as the molecule bends or as 
one CO bond is extended. The two intersections, implied in Fig. WCh), are not independent: 
They belong to a single 'CI seam' in the {Ri, R2) plane of the CO bond distances in linear 
OCO. The CI seam forms a closed loop as discussed in Sect. HI C of paper I and in Ref. 



In fact, the closed CI seam is a line of fivefold degeneracy which passes through the FC 
region. Calculations discussed below demonstrate that the closed seam substantially affects 
vibronic eigenstates located above the CIs in the strong absorption band at 133 nm. 
(3) The Herzberg- Teller (HT) effeciP' is seen as a strong coordinate dependence of the 
transition dipole moments (TDMs) with X. The TDM vectors in the principal axes system 
are shown in Fig. [3] as functions the CO bond length and the OCO angle.'^Zl Although the 
HT interaction neither requires nor involves a direct state intersection, its effect upon the 
intensity distribution in the spectrum is virtually indistinguishable from that of a CI in 
the FC region. Optical transitions from X^S+ to the states 1^S~, 1^11^, and l^A^ are 
forbidden in the -Dqo/i symmetry group. The TDM along the molecular axis, fiy, transforms 
according to the irrep F'^^ = a^, while the components fix and //^ are of F'^^" = F^"^ = iiu 
symmetry. Direct productP^El! Yf x F^ x F^^ of these irreps with the electronic symmetries 
Ff = S^ and F'^' = S~, II^, or A„ are never totally symmetric and all TDMs at the 
FC point vanish. Transitions become vibronically allowed in distorted molecule, and the 
TDMs directly depend on displacements in non-totally symmetric coordinates. The observed 
absorption between 120 nm and 160 nm is entirely due to this coordinate dependence. For 
linear deviations from Dooh, active in the transitions are the antisymmetric stretch i?_ = 
{Ri — R'i)l^2 (irrep cr+) and the bend aoco (irrep 7r„). In linear CO2, the components 
[ix,z are antisymmetric with respect to i?_ as illustrated in Fig. [3ta,c). The component 
liy grows linearly with aoco; while [ix,z are almost independent of the OCO angle [Fig. 
3[d,e,f)]. The TDM components corresponding to states 1^S~ and 1^A„ are easily identified 
at linearity [Fig. p[a,c)] by their small absolute values; the transition to l^II^ dominates 
in Cooj;. The coordinate dependence of the TDMs, violating both the Born-Oppenheimer 
and the FC approximat ions ,121251 jg commonly ascribed to the 'intensity borrowing' from 
an adjacent bright lender state. The nearest bright state which might contribute to the 
vibronically induced TDMs in Fig. p^ is the Rydberg state l^S^ lying almost 2eV above 
l^IIg. As deviations from -Doo/i grow, higher than linear expansion terms contribute, all TDM 
components become nonzero, and all states can be excited in the bent molecule [Fig.[3Fd,e,f)]. 
The A states are excited via the in-plane components [ly^z- The out-of-plane component [ix 
mediates excitation of the A' states. The symmetry of \ix,y,z in Fig. [3] implies that parallel 
(perpendicular) transitions excite vibrational states which are symmetric (antisymmetric) 
with respect to interchange, R\ -H- R^^ of the two CO bonds. 



In this paper, the absorption spectrum of CO2 is calculated using ab initio PESs of the 
first five excited electronic states, vibronic interstate couplings, and multistate quantum 
dynamics for the nuclei. The study was designed with three goals in mind: To reproduce 
the experimental spectrum, to assign the calculated absorption bands, and to examine the 
relative contributions of the above non-adiabatic interactions. These tasks are connected, 
because photoexcitation brings CO2 into the region of multiple degeneracies (Fig. |2]), and 
all low lying states, populated via the HT effect, might be interlinked through an intricate 
network of couplings created at the RT intersections and the CIs. A realistic assignment 
should reflect the extent of this mixing. The outline of the paper is as follows: The multistate 
quantum mechanical calculations of the absorption spectrum and resonances in dissociation 
continua of the five excited electronic states are summarized in Sect. [TTl The calculated 



spectrum is compared with experiment and assigned in Sect. Ill Individual rotoelectronic 



transitions are isolated in Sect. Ill A allowing 'coarse grained' electronic assignment of the 



two calculated bands and a qualitative assessment of the role of RT coupling. Section III B 
assigns vibronic progressions in the stronger high energy band and shows that nodal lines 
of resonance wave functions follow the topography of the upper adiabatic PESs near the 
closed CI seam. Vibronic assignments of the low energy band, much less definitive because 



of strong anharmonic and vibronic couplings, are discussed in Sect. IIIC Weak absorption 



of the state 1^S„ is the subject of Sect. HID Summary is given in Sect. IV 



II. QUANTUM MECHANICAL CALCULATIONS 

A. Molecular Hamiltonian 

The starting point for the calculations is the 5x5 molecular Hamiltonian written in the 
basis of the electronic states |A'), |n'), |S"), |n"), and |A"), diabatized with respect to CIs 
at bent geometries and in the PC region as described in Sect. IVB of paper I: 

Ho = T + S^V^S = T + V^ (1) 

with the nuclear kinetic energy matrix T and the potential matrix V^, related to the diagonal 
matrix of ab initio energies V^ via orthognal adiabatic-to-diabatic (ADT) transformation S 



and consisting of a 2 x 2 block of A' and a 3 x 3 block of A" states: 



Va' Vu'A' 
Vri'A' Vw 



\ 



V = Vs" Vs"n" Vj^'iA" ■ \^) 

\ Vs"A" yn"A" Va" / 

All matrix elements are smooth functions of molecular coordinates; couplings of the ac- 
cidentally degenerate states |S") and |A") to |n") are set equal, and the RT-like S"/A" 
interaction is neglected. This representation will hereafter be referred to as Cl-diabatic; by 
construction, the projection of L on the molecular figure axis, L^, in the Cl-diabatic basis 
has a definite value of A = 0, 1, or 4. The kinetic energy operator T is specified in the 
body fixed (BF) frame in Jacobi coordinates {R, r, 7) in the form, originally suggested by 
PetrongolcP^ and modified for numerical implementations by Goldfield et al.P^ 



T 



fR + fr + {Br + Br) % + 



Jl + Ll 



sin^7 



Bui J'- 2J 



i{Br + 5,; 



JzLz 

sin^7 



— 2BrJ:zLz 



+ 



BrJ+3^ + 5r J_Jh 



(3) 



Jacobi coordinates in this expression comprise the distance R between one oxygen atom 
and the center of mass of CO, the CO distance r, and the angle 7 between the vectors R 
and r. Operators Tr, T^, and T^ are the standard^ kinetic energies for these coordinates; 
the rotational constants Br = {2^rR^) and B,. = (2^.,r^) depend on the corresponding 
reduced masses ^iR^r- J denotes the total angular momentum, J = iVco2 + L (iVco2 is 
the angular momentum of nuclei). J^ and J± denote the BF z projection and the ladder 
operators of J, respectively. The BF z axis is assumed to run parallel to R ('R-embedding', 
Ref. [371); finally, j± = iJz — Lz) cot 7 =1= d/d'y. In the basis of Cl-diabatic states, the first 
bracket in Eq. ^ applies to the vibronic dynamics in either A' or A" symmetry block; the 
second bracket is the RT interaction between A' and A" states sharing the same A^; the 
thrid bracket is the Coriolis interaction which includes both intra-symmetry (through e.g. 
J+Jz terms) and inter-symmetry (through e.g. —J+Lz terms) couplings. 

Divergence at 7 — )■ 0, tt of the kinematic RT coupling matrix elements for states with 
A^ 7^ is removed by switching from the definite symmetry states |A') and |A") to the 



RT-diabatic eigenstates |A ) of the Lz operator at linearity,'^ 

|A^) = ^(|A')±.|A")); (4) 

Applied to Eq. (|2|, this transformation gives a full CI/RT-diabatic representation, with 
all kinematic couplings replaced with potential matrix elements. The matrix V^ loses its 
block diagonal form and comprises three groups of elements: (1) Diagonal CI/RT-diabatic 
'potentials' |(Va/ + Va")? (2) off-diagonal vibronic couplings due to CIs involving states with 
Ai 7^ A2, ^(Va'a' + Va"a"), and (3) off-diagonal rotovibronic couplings ^(Va./ — Va") due to 
RT effect. The state |S"), vibronic couplings for which are set equal to those of |A") and 
the RT effect is absent, becomes decoupled from other states in this approximation. 

The rotoelectronic basis in which the molecular Hamiltonian is set consists of direct 
products {|0j{fn) |A^)} of the electronic CI/RT-diabatic functions and the rotational basis 
functions IOa/q) expressed via Wigner D matrices; M denotes the space-fixed and Vt the 
body-fixed projection of J . In the actual implementation, this basis is further transformed 
by taking parity adapted linear combinationgSSESl Qf basis functions with ±|/| = ±(i7 + A) 
or ±|/| = ±(fi — A). Further, the Jacobi angle is discretized and the angular differential 
operators are set in the Gauss-associated Legendre discrete variable representation (DVR). 
Finite basis representation is kept only in the quantum numbers Vt and A. 

B. The initial state and the absorption spectrum 

CO2 in the rotational state (Jj,f2j) in X^S+ is described by the wave function 

vl/x(q|J.,a)=^o^''^^(q)|e5(^^J. (5) 

Here ^/'g" '(q) is a given vibrational state and q denotes three internal (e.g. Jacobi) coordi- 
nates. Since X is a S state, the initial total angular momentum Jj coincides with the nuclear 
angular momentum Nqo2- Centrifugal sudden (CS) Coriolis-free approximation is used in 
Eq. ([5| treating CO2 as a symmetric top with conserved helicity fij. Only excitations from 
the vibrational ground state are considered, which for linear CO2 implies Vti = 0, too. 

The state of CO2 immediately after photoexcitation is first set in the Cl-diabatic repre- 
sentation. TDMs with X are obtained from ab initio TDMs using the ADT matrix S: 

M = V^A'5 A^n'' /^S"' Mn"' Ma"] = ^ /^ =^ {l^2A'yl^3A'yl^lA"yl^2A"jl^3A") ' ("J 
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The resulting TDM 'vector' has five electronic components, each featuring up to two non- 
vanishing spatial projections in the principal axes frame. Possible rotational transitions in 
the optical excitation are Ji -^ Jf = Ji, Jj ± 1 and Qi = ^ Qj = 0,1 (parallel case) or 
fij = — 7- fi/ = 1 (perpendicular case). Parallel and perpendicular transitions are calculated 
separately. The electronic component of the initial excitation in the Cl-diabatic state A with 
the angular momentum Jf and the BF projection Qf reads as 

Hq\jf,nf)A = {16, mjf^j) (4)5<''°(q) \QZ+M)n,) , (7) 

Here (■,■ | ■) is the Clebsh-Gordon coefficient, 6 = Qj — Qi = or 1, and (/iA)^ are the 
spherical projections of the TDM vectors ^^a: {ij.)s=o = l^z and {iJi)s=+i = — (fJ^y — ij^x) / V^- 
This expression is valid for both parallel and perpendicular transitions. The only difference 
between the two is that the initial state in a parallel transition is real, while in a perpendicular 
transition it is imaginary. For each type of transition, the total initial wavefunction $(q| Jj) 
is a linear combination of all $(q| Jj, Vtf)\ with the allowed values of A and {Jf, Vtj)^ Initial 
state with each Jf is propagated separately. The corresponding cross section, averaged over 
the initial space fixed projections M , is proportional to the expectation value of the Green's 
function of the time-independent Schrodinger equation,'^ 

a{E,^\Jf ^ J,) = ^($(q| J;)|ImG+(EpO|$(q| J/)) , (8) 

where the photon energy E^Yl is measured from the zero point energy level EQ{Ji,Qi = 0) in 
X. The total cross section is the sum over all allowed J/: 

a{Eph\Ji) = V a{Epipf ^ Ji) . (9) 

The calculations including all Qf blocks for each Jf and thus fully including Coriolis coupling 
between adjacent helicity blocks were performed for Jj < 2. 

In order to compare the calculation with experiment at T = 190 K, a thermal absorption 
spectrum has been constructred by averaging a{Ep^\Ji) over the normalized Boltzmann 
population Pj^T) of the initial states Ji, 

a{Ep^,T) = J2PjXTME,^\Ji)- (10) 

Ji 

Because of large Ji < 45 involved in these calculations, CS approximation was applied, 
the Coriolis coupling neglected, and each initial state $(q|J/,fi/) processed separately. For 
Ji < 2, the spectra obtained with and without CS approximation are in excellent agreement. 
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In practice, calculations of the absorption spectrum via Eq. (pj) are performed in several 
steps. The Cl-diabatic initial state $(q| J/) is transformed, first to the CI/RT-diabatic and 
then to the parity adapted representation. The last transformation is dropped in the CS 
approximation. Next, the outgoing boundary conditions are approximately imposed in the 
asymptotic region by augmenting Hq with a complex absorbing potential —iW (Re W > 0)*^ 

H = Ho-ziyi. (11) 

Here 1 denotes the 5x5 unit matrix, and one coordinate function W is used in all electronic 
channels. In the last step, the vector lm.G^{E) $(q| Jj) is approximated by the expansion 

Alitor 

ImG+(E)$(q|J;) ~ J][Re6„(E)]$„, (12) 

n=0 

in which &„(£') are the usual energy-dependent coefficients,Sil while the vectors $„ are found 
from a modified Chebyshev recursion relation due to Mandelshtam and Taylor.l^ Chebyshev 
autocorrelation coefficients c„ = ($o|*&n) with the initial state $0 = "^(qI^/) are calculated, 
and the absorption spectrum is reconstructed from {c„} using the Eq. ([s]). The number of 
expansion terms, iViter, is one of the convergence parameters of the calculation. 

In order to assign diffuse peaks in the spectrum, filter diagonalizationpES has been per- 



formed. Eq. (12) is used to generate a set of basis vectors {ipk}, 

^fc = ImG+(Efc)$o, (13) 

localized around nodes {Ek} of an equidistant energy grid covering a given window 
[-Emin, -f'max] (500 cm^^ widc wiudows were used in most calculations). The Hamiltonian 
matrix (ipklHlipk/) is set in the orthogonalized basis and diagonalized. The resulting com- 



plex eigenstates are resonance states of the Hamiltonian (11) with energy En and width F^. 



The eigenfunctions \l/„ are five component vectors known in any representation used in the 
calculations. For example, the Cl-diabatic eigenstate, ^„ = (\I^^ , \I'^ , ^^ 5 ^n 5 ^n ) 1 
can be rewritten in the Cl-adiabatic representation using the ADT matrix via S"^^„. 

C. Numerical details 

Quantum mechanical calculations have been performed with several sets of the diaba- 
tized PESs constructed using slightly different fitting procedures for ADT parameters as 
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described in Sect. IVB of paper I. All diabatic sets give indistinguishable absorption spec- 
tra. Hamiltonian is set in DVR"^^ with the grid comprising 150 potential-optimizecP^ points 
in i? G [1.5, 7.0] oq, 100 potential-optimized points in r G [1.5, 6.5] Oq, and N^ = 100 Gauss- 
associated Legendre quadrature points in angle. The grid is contracted by retaining only the 
points with potential energy below V^ut = 12.5 eV above equilibrium in X. Niter = 15 000 
Chebyshev iterations are sufficient to converge the absorption spectrum. The complex ab- 
sorbing potential —i W is set via a coordinate dependent damping function 7,E21 which is 
zero in the inner parts of the five adiabatic potentials and grows quadratically as a function 
of CO- ■ ■ O distance in the asymptotic regions. The quadratic polynomial 7((5) {Q = R 
or r) is described in Refs. [39landB6l It becomes non-zero 1.1 oq before the grid edge and 
reaches the strength of 0.15 at the edge. 

Temperature dependent absorption cross section is based on spectra calculated for all 
initial Jj between and 45, the range sufficient for temperatures up to at least T ~ 250 K. 

Filter diagonalization is performed in several overlapping energy windows covering the 
photon energy interval from 60 000 cm~^ to 80 000 cm"^. Calculated resonances are broad, 
with the average width of r„ > 50 cm~^, and require 25 000 Chebyshev interations for 
convergence. As discussed in Sect. |IIIB the calculated resonances allow a qualitatively 



accurate reconstruction of the absorption spectrum. All quantum mechanical calculations 
in this work are performed using the program package 'Poly Wave '.^3 

III. RESULTS 

The absorption spectrum, thermally averaged at T = 200 K, is compared with the exper- 
imental cross section in Fig. [lib). The calculated spectrum has a characteristic two-band 
shape, and the band centers at 8.43 eV (147.2 nm) and 9.28 eV (133.6 nm) are accurate within 
0.05 eV. The maximum cross sections in the two bands, 5.0 x 10^^^ cm^ and 9.5 x 10^^^ cm^, 
are about 15% below the experimental values, but the observed 1:2 intensity ratio between 
the weaker and the stronger band is accurately reproduced. In both spectra, structures in 
the weaker 148 nm band are substantailly more diffuse and less regular than in the bet- 
ter resolved stronger 133 nm band. Large discrepancies are found in the spectral widths 
(FWHMs). For the low and the high energy band, FWHMs are 7500 cm~^ and 6100 cm~^, 
i.e. off by 6% and 18% as compared to the respective values of 7000 cm~^ and 7400 cm~^ 
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obtained from deconvolution of the experimental spectrum of Ref. [161 The exaggerated gap 
between the two calculated bands is one of the reasons for this deviation. The gap depth, 
only weakly dependent on diabatization details,'^ is sensitive to the TDM derivatives with 
respect to molecular coordinates at the FC point. In several test calculations, these deriva- 
tives — the manifestations of the HT effect — were artificially modified, and the angular 
derivative d fiy/daoco had the largest impact on the band intensities and the gap size. 

The impact of thermal averaging can be assessed by comparing the 200 K spectrum in 
Fig. [IVb) with the 'OK' spectrum, a{Epii\Ji = 0), calculated for a single initial rotational 
state Ji = and shown in Fig.llFa). The differences in positions and maximum intensties are 
marginal: The 200 K spectrum is red shifted by ~ 100 cm^^ and its intensity is about 12% 
higher. Boltzmann averaging primarly affects the extent of spectral fluctutations which in 
the 200 K spectrum is substantially attenuated. The reason is the inhomogeneous broadening 
in the spectrum 'assembled' from an incoherent average over many spectra for different 
rotational states Jj. Narrowly spaced and erratically fluctuating lines in the low energy 
band are more amenable to thermal 'washing out' than the well resolved peaks in the high 
energy band, and the temperature effect is most conspicuous below 9 eV. 

A. Electronic assignments 

The currently prevailing assignment, attributing the 148 nm band to 1^A„ state and 
the 133 nm band to l^H^ state,'^^'^^'^ is somewhat ambigous because these states cross in 
the FC region. Moreover, this assignment clearly refers to linear geometries because CI- 
diabatic picture is explicitly selected as a reference. Quantum mechanical approach offers 
direct means to clarify the electronic origin of each band. For example, the assignments in 
Ref. [T] were determined using population patterns of the electronic components of resonance 
wave functions ^„. The results indicate that it is the Cl-adiabatic representation which is 
adequate for the assignment of resonances, implying a strong vibronic coupling scenarioP^ 
The 2^A' and l^A" adiabatic components are strongly populated below the CIs, while the 
components 3^A' and 3^A" become dominant above the CIs. However, the population analysis 
works best if the electronic states involved are all bound or all unbound. If, as in the present 
case, the molecule is excited into the energy range in which the linear states 3M'/3M" are 
effectively bound, while the bent states 2^A' /l^A" are dissociative, the results might become 
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sensitive to the precise definition of the 'inner region' used to evaluate electronic populations. 
The Cl-adiabatic assignments of Ref. [T] can be corroborated using short-time dynamics 
— a natural approach to establish the electronic labels which refer to the largest energy 
scale (i.e. the smallest time scale) in the molecule. To this end, two auxiliary diabatic TDM 
vectors are constructed using the ADT of Eq. (|6]), one exclusively from the TDMs of the 
bent adiabatic states 2M'/lM" (termed ^^21) ^^^ ^^e other exclusively from the TDMs of 
the linear adiabatic states 3M'/3M" (termed ^^33): 

M21 = S^(m2A', ,MiA",0, f (14) 

M^3 = S^( ,M3A', ,0,M^^„)'^- (15) 

For each auxiliary diabatic TDM, the initial excitation, constructed using Eq. ([t]) for Jj = 0, 
is propagated under the full Hamiltonian for 500 Chebyshev iterations corresponding to a 
time interval of ~ 0.1 fs. The resulting low resolution spectra are shown in the upper part 
of Fig. liFb). They provide a pictorial decomposition of the full spectrum a{Ep\^\Ji = 0) 
into the components due to excitations via ^21 ^^'^ Mss- Although not quantitative, the 
result is clearly consistent with the previous assignment of the 148 nm band to the bent pair 
2M71M" and the 133 nm band to the linear pair 3^A'/3^A". 

The Cl-adiabatic assignment is also motivated by the shape of the initial excitation along 
the bending angle, sketched in Fig. ^a) together with the cuts through the calculated PESs. 
The dashed red line is the ground vibrational state in X. Its maximum is shifted by 4° away 
from linearity due to the zero-point bending vibration plus the weighting imposed by the 
volume element. Thus, even before the TDM with the excited states is applied to ^Pq' (q), 
the FC geometry corresponds to aoco = 176°, and the gap between the adiabatic linear 
and bent excited states is ~ 0.6 eV. The solid red line in Fig. Wia) exemplifies the influence 
of the TDM. Shown is the initial excitation created in the state 3^A' via the TDM fiy (the 
initial excitation in 2^A' looks very similar), fiy strongly depends on angle because of the 
HT effect, as shown in Fig. Isle). As a result, the non- vertical transition bends the molecule 
even stronger, shifting the probability maximum to aoco = 173°. The linear /bent gap at 
this angle becomes 0.86 eV and matches the energy separation between the maxima of the 
two absorption bands (0.85 eV). In fact, the adiabatic gap is a primary control handle of the 
positions of the band maxima: If the original potentials are modified to reduce the gap size 
step by step, the bands gradually approach each other and eventually merge. 
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While vibronic couplings and the HT effect directly shape the UV absorption of CO2, 
the role of the rotoelectronic RT coupling is more subtle. The lower part of Fig. |4Fb) shows 
the parallel, cr||, and perpendicular, a±, components of the full spectrum cr(i?ph| Jj = 0). In 
the parallel transition, both Qf = and Qf = 1 states are initially populated in the A' 
symmetry states, but 86% of the population is in the Qf = block reached via the strongly 
angular dependent TDM fiy. RT coupling is zero for Qf = 0; it is active only in the almost 
empty flf = 1 block, reached via fi^, and is promoted by the Coriolis-induced population 
transfer from Qf = OtoQf = l. If the molecule dissociates exclusively out oi Qf = 1 states, 
Coriolis transitions control the dissociation rate. Dissociation of HCO in the A^A" state is 
a well known example.'^^'^ In contrast, dissociation of CO2 is fast even for Vlf = and the 
Coriolis interaction is inefficient on the characteristic time scale of rdiss ~ fi/Tn ~ 100 fs. 
This is also the reason why CS approximation has little effect on the parallel spectrum a\\. 

The perpendicular transition terminates in A" states and is realized via TDM /ij,. The 
spectrum a± in Fig. |4Fb) consists of two widely separated almost equally intense bands 
featuring apparently regular progressions of diffuse peaks. The reason for the distinct reg- 



ularity of the 148 nm band in a± will be discussed in Sect. IIIC The rotational state 
(J/ = l,fi/ = 1) is exclusively populated, and the RT coupling is operative in the per- 
pendicular transition. However, a± is substantially smaller than a\\ for all but the highest 
photon energies: The TDM fj.^., which is similar to /x^, depends weakly on aoco (see Fig. |3]), 
and it is excitations from bent geometries which dominate the spectrum. As a result, the 
RT coupling, active in a±, has little influence on the sum a±_ + (J\\- 

The state ^S^, decoupled from other states in this calculation, is excited in a perpen- 



dicular transition via fix- Its absorption spectrum, discussed in Sect. |IIID[ consists of a 
series of sharp narrow lines, with the peak intensity reached in the gap between the two 
main bands. The TDM with X is extremely small in the FC region, and the absorption 
cross section does not exceed 5 x 10~^°cm^. As a result, the contribution of the ^S~ state 
is hardly discernible in the full spectrum. Test calculations register minor intensity gain if 
weak vibronic coupling to other states is included. 
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B. Vibronic assignments of the high energy band 

The intense lines in the experimental 133 nm band [Fig. [iFb)] are spaced, on average, 
by ~ 630 cm"^. This major progression is shown by brown sticks in Fig. p[a). Another 
progression can be recognized as minor peaks or shoulders above 72 000cm~^. The mi- 
nor peaks, also shown as sticks in Fig. p[a), are red shifted by ~ 250 cm~^ with respect 
to the major ones. Above 77000cm^^, the diffuse bands consist of two clearly separated 
approximately equally intense peaks. Despite its simple appearance, even the major pro- 
gression resists straightforward assignment because the spacings between the adjacent peaks 
are not constant. Instead, they obey an intricate rise- and- fall pattern in which intervals of 
frequency growth (for example, 73 000cm~^ — 75 000cm~^ or 76 000cm^^ — 78 000cm~^) are 
interrupted by dips (for example, near 73 000 cm~^, 75 300 cm~^, or 79 000 cm~^). Rablais et 
al. suggested that the progression is due to bending excitations in the near-linear molecule, 
with the perturbations originating from the Coriolis interaction between the vibrational and 
rotational angular momenta.'^ The dynamics calculations suggest that the major progression 
is built on synchronized excitations of the two CO bonds. 

In the calculated spectrum [Fig. fltb) andllFa)], the major and minor progressions are 
qualitatively similar to the experimental ones, but numerically different. For example, the 
spacings between major peaks follow the same rise-and-fall pattern, but their average is 
merely ~ 580 cm~^. Next, the minor progression, also seen in the calculations, is blue- 
shifted by ~ 240 cm~^ against the major one. Around 76 000cm~^, the intensity of the 
minor peaks grows substantially and the calculated spectrum shows several strong lines giv- 
ing it a congested appearance. Above 77000 cm~^, most calculated peaks are split into two 
components in agreement with experiment. Thus, the achieved accuracy is not spectro- 
scopic, but the calculated spectrum represents a reasonable starting point for the analysis 
of molecular motions behind the diffuse bands. The assignment, presented below, is based 
on the properties of metastable resonance states calculated specifically for the parallel com- 
ponent cr|| of the full spectrum. The principal progressions are similar in a{Ep^\Ji = 0) and 
a\\{Eph\Ji = 0) [cf. Fig. E[a) and (b)]. Differences visible above 78 000cm~^ are addressed 
at the end of this section. 

Positions, En, of resonance states excited in a parallel transition {Ji = 0, fi^ = 0) — ?■ {Jf = 
l,fij = 0), calculated using filter diagonalization in the CS approximation, are illustrated 
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with sticks in Fig. W[c). Each resonance is further characterized by its width r„ and the 
intensity /„ = |($(Ji, fij)|^„)| . The resonance absorption spectrum, 

.,.(EpHiJ,.n,)=x:/,.^^— 0^^. (16) 

is compared with o"|| in Fig. W[c). Resonances are expected to accurately describe absorption 
spectra consisting of narrow isolated lines.'^ However, despite a substantial fast direct contri- 
bution, the incoherent sum of Lorentzians in o"res is in qualitative and — above 73 000cm~^ 
— even in quantitative agreement with (T|| , both in terms of positions of the diffuse lines and 
the intensity of the background. 

Nodal patterns of resonance wave functions provide vibronic labels for the diffuse peaks. 
All resonances with r„ < 150 cm~^ were screened, but the discussion will be limited to states 
with the largest intensities. Their positions are shown with sticks in Fig. |4Fd), and their 
widths range from 30cm~^ to 80cm~^. Relevant for the assignment is the 3^4' electronic 
component of the ^„ vector in the Cl-adiabatic representation. Indeed, the 133 nm band 
is due to the linear adiabatic states 3^A'/3^A". Moreover, the considered parallel transition 
populates only A' states, with the population transfer into the A" block suppressed for 
Qf = without Coriolis coupling. Direct dissociation in the state 3^A' is hindered by 
a ~ leV high barrier [see Fig. 1 of paper I], and the 3 A' component of ^„ describes a 
vibrational level trapped in the upper adiabatic sheet of the cone and non-adiabatically 
coupled to the dissociation continuum of the lower 2 A' sheet (cf. Ref. [52]) . 

States in the major progression, which can be followed from the origin near 72 700cm~^ 
up to at least 78 000 cm~^, are illustrated in the left column of Fig. [sj The 3D wave functions 
are viewed in the (i?i, R2) plane against the contour map of the 3 A' potential containing the 
closed CI seam. Their assignment in terms of three quantum numbers (f^, Vr, Vb) is indicated 
above each frame and the respective combs are shown in Fig. |4Fd). These quantum numbers 
reflect three features shared by the states in the main progression. First, they have no 
nodes along the bending direction except at aoco = 180°; most probable angles lie between 
173° and 177° and the states are essentially 'flat', i.e. localized in the (-Ri,i?2) plane. Such 
resonances are labelled by Vf, = 0. Second, these states have no nodes along the direction 
normal to the CI seam, and their 'radial' quantum number is zero, Vr = 0. Third, these 
states decribe consecutive excitations around the CI loop. In each successive state, two new 
nodes are 'stringed' on the CI seam 'wire', and the quantum number for this mode, v^, 
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counts nodes along one side of the seam. All states in Fig. [5] are symmetric with respect to 
interchange of -Ri and R2, as expected for parallel transitions; antisymmetric states, with a 
node along the Ri = R2 line, remain 'invisible'. The potential along the seam loop is rather 
smooth, and the average frequency in the Vcf, progression is w^ ~ 580 cm~^. However, energies 
in the 'angular' v^ progression are not equidistant (see below). The quantum numbers Vr and 
Vb correspond to the tuning and the coupling modes of the CI, respectively. The adiabatic 
potentials along these directions are cusped, and the vibrational frequencies Ur and uji, are 
expected to be much larger than u^. 

A simple two-dimensional (2D) reference model, which illustrates the impact of the closed 
CI seam on the eigenstates and the energy spectrum of the full Hamiltonian, consists of a 
single adiabatic 3M' PES in the (/2i,i?2) plane with fixed aoco = 179°. The 2D spectrum 
has been calculated using filter diagonalization, and examples of the 2D eigenstates are 
shown in the middle column of Fig. Isj They can be assigned two quantum numbers {v^^ f t-)2D 
measuring excitations along the angular seam mode and the radial stretching mode r. As 
in the full Hamiltonian, the states have either a maximum or a node along the line Ri = R2- 
Symmetric wave functions in the progression (f(/,,0)2D in Fig- [S] can be directly compared 
with states in the major progression (t;(^,0,0) in the left column. The 2D states are tightly 
localized around the CI loop, emphasizing that the quantum number v^ describes hindered 
rotation in the (i?i, -R2) plane. Similarities in the wave functions are strong enough to make 
the 2D model useful in disentangling complicated nodal patterns in the full calculations. A 
good case in point is the progression (f^, 1)20 built on one quantum of radial excitation and 
shown in the lower panels of Fig. [6J The frequency of the radial mode, cOr = 2150 cm~^, is 
almost four times the hindered rotation frequency. Resonance anharmonic couplings mix the 
states {V(j,, 0)20 and {v^ — 4, 1)20, but their wave functions remain clearly recognizable in 2D 
(compare states marked with stars in Fig. Island |6]). The progression of states with f ^ = 1 is 
found in the full calculation, too (Fig. [61 upper panels). The fundamental radial frequency, 
Ur = 1560 cm~^, is lower than in 2D. Nevertheless, the anharmonicity is strong and many 
states {v^, 0, 0) and (f ^ — 4, 1, 0) are mixed making their assignment difficult. In fact, the 
states in the upper panels of Fig. |6] were identified using their 2D analogues. Moreover, 
most states {v^, 1,0) are discenrible merely as satellites to the major peaks {v^,0,0); their 
intensity increases only above 78 000cm~^ [see Fig. p[d)]. As a result, many members of 
this progression are missing, but the gaps can be filled in by the 2D model (cf. Fig. |6]). 



17 

The third progression clearly visible in the calculated spectrum consists of resonances 
with a node at a bending angle aoco < 180°. Such states are assigned {v^,Vr, 1). The state 
(0,0, 1) lies at 75 300cm~^, so that the bending frequency is ut, = 2290 cm~^ The intensity 
in this progression quickly grows with energy. Of the three strong bands, 'congesting' the 
spectrum between 75 000cm~^ and 76 000cm~^, two are due to states (2,0, 1) and (3,0, 1); 
only one state, (6,0,0), has no bending excitation. Above 76 000cm~^ the states (f0,0, 1) 
have the largest intensities. Their wave functions, shown in the right column of Fig. |5| 
are similar to the respective states with vi, = and can again be analyzed in the (i?i,i?2) 
plane. Bending excitation introduces an additional 'layer' of nodes and enhances mixing 
between the combination states {v^,0, 1) and (f^ — 4, 1, 1) and both assignments are given 
in some panels. Above 78 000cm~^, the assignments become even less unique. Several 
factors contribute: Resonances substantially broaden, states with Vb > 1 are found, and 
mixing among the three vibrtaional motions increases. Conservative assignment associates 
the simple looking high energy bands with progressions {v^, 0, l)/(f0 — 4, 1, 1) and {v^, 1, 0), 
although this selection is by no means unique. 

The shape of the lower adiabatic potential in the plane of the tuning and the coupling 
modes of a CI is often referred to as a Mexican hat, especially in the context of the E x e 
Jahn- Teller effect.^ A closed CI seam discussed in this paper adds another hat to the rack: 
A 'Mongolian hat', sketched in Fig. [5} is a piece of headwear which describes best the shape 
of the upper adiabatic PES in the {Ri,R2) plane containing the CI seam. The hindered 
rotation excitations are winding along the hat brim (the CI seam), while consecutive radial 
excitations extend towards (and eventually engulf) the top of the hat. 

The closed CI seam leaves a clear mark in the calculated absorption spectrum. Frequen- 
cies in the progression of pure f ^ excitations are shown in Fig. [tI in panel (b) for the full 
calculation and in panel (c) for the 2D model. Consider the model first. Energy intervals 
grow between Ep^ = 73 000cm~^ and 76 000cm~^, and this is consistent with a pseudoro- 
tational motion around the seam and the associated quantization with E^ ~ v'i. At the 
same time, the frequency in the (f<^,0)2D progression is not monotonic, and there are two 
dips at low and high energies. Minimum in a progression frequency is often an indication 
that a critical energy in the potential (a barrier top or a saddle)^ or in the phase space 
(a separatrix) has been reached.'^^'^ Possible sources of the rise-and-fall pattern in the 2D 
model are illustrated in the lower panels of Fig. [7] showing one-dimensional cuts through 
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the PES of the adiabatic 3M' state, taken along the seam [panel (d)] and across the seam 
[panel (e)]. Angular coordinate (f) along the seam in (d) is chosen such that (f) = 180° for 
Ri = R2 = 2.24 ao] the 'radial' coordinate in (e) runs along the antisymmetric stretch R_, 
and the origin i?„ = is the point Ri = R2 = 2.5 ao- The potential along the seam has two 
global minima separated by a local one at = 180°. The first frequency dip (gray color) in 
the (f 0, 0)20 progression involves the excited state in the global well and the ground state 
in the local well. The second frequency dip (yellow color) involves the states located near 
two types of critical points: The barriers, separating the local and the global minima along 
the seam [(f) = 135° and 225° in panel (d)], and the top of the Mongolian hat, separating the 
global minima across the seam [i?_ = in panel (e)]. 

The principal topographic features, captured by the 2D model, shape the spectrum in the 
full calcultion, too. The frequency in progression {v^, 0, 0) in Fig. [Wb), although numerically 
different from the 2D case, grows steadily up to -Eph ~ 75 000cm~^ reflecting consecutive 
pseudorotational excitations along the CI seam. The frequency dip at high energies indicates 
that the top of the Mongolian hat has been reached, while the low energy dip is absent 
because the origin of the progression lies above the bottom of the local minimum at = 180°. 

The accuracy of the ab initio PESs does not allow a line by line comparison with the 
measured spectrum. To which extent the above assignment is applicable to the experiment? 
Energy intervals in the major progression measured below Ep^ = 78 000 cm~^ are shown in 
Fig. ^a), and the rise-and-fall pattern, inherent to the calculated spectra, appears to be 
distinctly recognizable, including a stretch of pseudorotational excitations and dips impying 
potential saddles. This observation should not be overinterpreted — shown is one progres- 
sion, and this progression is rather short. If, however, the agreement with the calculations is 
more than a mere coincidence, the 133 nm band would be the first documented spectroscopic 
observation of a closed CI seam. 

I would like to conclude the discussion of the 133 nm band by pointing out an alternative 
explanation for the double peaked lines above 78 000cm~^. Their assignment in terms of 
progressions (f(^,0, 1) and (f,^, 1,0) has been given only for the parallel component of the 
total spectrum. Decomposition into o"|j and a± at the bottom of Fig. W[h) demonstrates that 
intensities in the two transitions become comparable at high photon energies. In fact, each 
double peaked line in the total spectrum consists of a low energy parallel peak (green line) 
overlapped with a high energy perpendicular peak (red line). 
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C. Vibronic assignments of the lo^v energy band 

Assignment of the 148 nm band is a challenging problem: Diffuse structures are weak, 
irregular, and overlapping. Substantial bending excitation has been suggested for this 
bancP'^ originating from states 2^A' and iM" with strongly bent equilibrium geometries.'^^'^ 
The present calculations demonstrate that actual molecular motion is difficult to define in 
terms of the familiar normal modes: Energy above equilibrium is high and the anharmonic 
and vibronic couplings are substantial. The results summarized in this section at best pro- 
vide an initial point from which a genuine analysis could be undertaken. 

In the experimental spectrum, almost every diffuse structure below 72 000 cm~^ comprises 
several transitions visible as a group of peaks or as a shoulder [Fig. [l]^b)]. Diffuse struc- 
tures in the calculated spectrum look similar, especially after thermal averaging, although 
irregularities in the energy intervals, widths, and intensities preclude a meaningful matching 
with experiment. In order to enable comparison, line subsets are selected in both spectra. 
In the experimental spectrum, centers of the diffuse peaks with the largest intensities are 
chosen as indicated in Fig. |4Fa). In the calculated spectrum, only the parallel transition 
(Jj = 0, fij = 0) — 7- (J/ = 1,^2/ = 0) is initially considered, and resonance states ^n with 
the largest intensities are taken to represent diffuse peaks in a^ [Fig. |4](d)]. 

Energy intervals between the selected peaks are shown in Fig. M (-Eph > 60 000cm~^). 
Experimental peaks are spaced by 400 cm~^ — 700 cm~^. Spacings between the intense res- 
onances vary in similar limits although the agreement is only qualitative. Moreover, visual 
assignment of the wave functions fails for the selected resonances. A typical example is 
given in Fig. |9] which depicts the 2A' adiabatic component of a state lying at E'ph = 8.38 eV, 
i.e. ~ 3eV above equilibrium of the 2A' state and ~ leV above dissociation threshold. 
The large excess energy is the reason for the bizzar nodal patterns of the 2A' resonance 
components, and — although wave functions of different states in Fig. |4](d) have different 
nodal patterns — their complexity is beyond rationalization in terms of progressions. 

An alternative 'bottom-up' approach in this case is to study the spectrum of the bent 2M' 
state from the ground state upwards and and to follow the characteritic frequencies with 
growing internal energy. Bent electronic states have C2v equilibria (see Table III of paper I), 
and the relevant vibrational quantum numbers are the symmetric stretch Vg, the bend Vb, and 
the antisymmetric stretch Va- AH wave functions in the full calculations, both resonance and 
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bound, were screened in order to find states, 2A' adiabatic components of which represent 
pure excitations in the three vibrational modes. Examples of pure stretching and bending 



wave functions are shown in Figs. 10 and 11, respectively. The progression {vs,0,0) is the 



simplest: Consecutive wave functions have parallel nodal lines (normal to which indicates 
the normal mode direction), and the symmetric stretch frequency Ug ~ 1400 cm~^ is mildly 
anharmonic and depends smoothly on energy (see Fig. ^ Ep^ < 60 000cm~^). In the pro- 
gression (0, 0, Va), only states with even Va are visible, because the transition is parallel and 
both ^y and /i^ are symmetric with respect to interchange of -Ri and R2. The antisymmetric 
stretch frequency, Ua ~ 900 cm^^, grows at low energies, passes through a maximum near 

notica- 



49 000cm~^ and then decreases (Fig. ^; nodal lines in the lower panels of Fig. 10 
bly bend in the {Ri,R2) plane with growing Va- The anharmonicity is the strongest in the 
bending progression (0,ff,,0). While the direction of the bending vibration is rather well 
preserved, its frequency in Fig. N varies from ujb{A') ^ 550 cm~^ down to 380 cm^^, experi- 
encing a slow rise followed by a rapid decline above 49 000cm~^. The mode coupling and 
the density of states grow substantially above i^ph = 52 000cm~^, rendering further search 
for pure excitations unfeasible. This leaves a 10 000 cm~^ wide energy gap between the last 
assigned state and the resolvable red edge of the 148 nm band. In view of the capricious 
anharmonicities, an extrapolation over this gap would seem unreliable for any mode but the 
symmetric stretch. However, if it is granted that the two lowest frequencies, Uf, ~ 350 cm~^ 
and Ua < 900 cm~^, keep their values up to and above 60 000cm~^, one would expect the 
diffuse peaks to reflect strongly excited bending and antisymmetric stretch motions. 

One major drawback of the above discussion is the missing experimental counterparts of 
the calculated vibrational frequencies: The measured intensities are low and the rotational 
structure is unresolved. Truly exceptional in this respect is the high resolution study of the 
wavelength range 172 nm — 198 nm, performed by Cossart-Magos and co-workers^ who were 
able to detect many absorption bands with rotational structure and to isolate two vibra- 
tional progressions between 51000cm~^ and 57000cm~^. This energy range overlaps with 
the energies in which pure excitations in the quantum mechanical calculations still could be 
detected. The progressions were assigned to the electronic state IM2, i.e. the bent state 
iM", excited in a perpendicular transition from X. Frequencies in the experimental pro- 
gressions are shown with open triangles and squares in Fig. |8j open circles are the frequency 
shifts between states belonging different progressions. In Ref. UHl both progressions were 
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interpreted as consecutive bending excitations differing by one quantum of the symmetric 
stretch, giving experimental estimates of u^ ^ 600 cm~^ and Ug ~ 1370 cm~^. Relevant 
theoretical progressions are shown in Fig. |8] with red symbols. They are obtained by cal- 
culating eigenstates in the perpendicular transition ( Jj = 0, Qi = 0) — )■ (J/ = l,Qf = 1) 
in many energy windows and looking for the wave functions ^„, 1^" adiabatic components 
of which represent consecutive bending excitations. Because fix{Ri,R2) = —fJ'x{R2,Ri), all 
excited states involve odd number of quanta in the antisymmetric stretching mode Va- The 
progression shown in the lower part of Fig. [s] (red squares) is (l,fb, 1), and the frequency 
shifts with respect to the progression (0, v^, 1) are given in the upper part of Fig. |8] (red 
circles). The calculation clearly confirms the experimental assignment. Agreement with 
the experiment is quantitative, especially for the bending frequency Ubi^A") in the (l,f6, 1) 
progression, which lies within 10 cm~^ of the experimental values. Somewhat less accurate is 
the symmetric stretch frequency shift and its energy variation. Most probably, this implies 
that the observed progressions are built on higher symmetric stretch excitations than Vg = 
and Vs = 1. The calculated intensity generally grows with Vg, but sufficiently long bending 
progressions with Vg >2 become tough to identify. 

The bending progressions in the parallel and perpendicular transitions shown in in Fig. 
Is] are strikingly different: ujb{A") in a± smoothly depends on energy and is only slightly 
anharmonic, while Ub{A') in o"|| has a strong non-monotonic energy dependence. This is 
surprising because the equilibria of the bent states 2^A' and iM" are similar. The sole 
difference is in the equilibrium angle, 127° in l^A" vs. 118° in 2^A', but it can hardly be 
made responsible for the observed discrepancy in the bending frequencies. Instead, the strong 
energy dependence of uJb{A') stems from another property of the 2^A' PES along the bending 
angle: In the calculation, it is this state to which the carbene-type bent OCO minimum at 
70° belongs after local diabatization of bent CIs. This minimum at aoco = 73° is is clearly 
seen in the contour maps in Fig. 11 Near E^i^ ^ 48 000cm~^, the ground vibrational 



state appears in this local OCO minimum shown in the left lower panel in Fig. [IIJ Above 
50 000 cm~^, the bending states in the local and in the global equilibria become connected via 
a saddle point, the main bending progression includes states partially delocalized between 



the two wells [right lower panel in Fig. 11 , the configuration space available to bending 
vibrations suddenly expands, and it is this energy range in which Ub{A') starts to decrease. 
In contrast, the l^A" state, diabatized at bent CIs in a different symmetry block, loses the 
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OCO bent minimuin which becomes reassigned to state 3M". As a consequence, the energy 
dependence of uJi,{A") is monotonic and inconspicuous. 

The difference in bending progressions in the states 2^A' and l^A" and its connection 
to the local OCO minimum along aoco allows one to explain the difference in the diffuse 
structures in the parallel and perpendicular excitations of the 148 nm band in Fig. |4](b). 
A nearly constant frequency uJb{A") correlates with the regular progression of peaks in a±. 
On the other hand, the bending frequency Ub{A') in the parallel excitation is small and 
strongly energy dependent — and the diffuse peaks in a\\ are narrowly spaced, irregular, and 
overlapping. Simultaneous analysis of both types of transitions indirectly confirms that the 
observed absorption lines are indeed due to strong bending excitations: Qualitatively dif- 
ferent bending potentials support qualitatively different spectra. Moreover, the 'congestion' 
in the 148 nm parallel band is interpreted as a direct consequence of the carbene-type OCO 
'cyclic' minimum in the bent 2M' state. 

It is worth reminding the reader that the question of electronic assignment of the carbene- 
type OCO goes beyond the diabatic states 2M', iM", or the ground electronic state in which 
it was discovered by Xantheas and Ruedenberg.!^ In paper I, the origin of this minimum 
was traced to the states A^A' and A^A" correlating at linearity with the Rydberg state 2^A„. 
The influence of the carbene-type OCO on the diffuse 148 nm band, uncovered in the above 
calculations, suggests that a systematic analysis of this band has to include the first four A' 
and A" states properly diabatized and fully interacting at multiple bent CIs. 

D. Vibrational assignments in the state ^S^ 



The spectrum of the state 2^A" (^S^ at linearity) is shown in Fig. 12 It stretches from 
65 000cm~^ to about 83 000cm~^ and consists of several series of narrow peaks. All vibra- 
tional states, excited in this transition, lie below 92 430cm~^, i.e. the threshold energy of 
the channel 0(^P) + CO(a^n) with which the state ^S~ correlates. Predissociation is purely 
non-adiabatic and therefore inhibited in these calculations. Because of the extremely small 
TDM Hx, the intensity of this band, excited in a perpendicular transition from X, does 
not exceed 6 x 10~^°cm~^. The frequencies of the bending, symmetric and antisymmetric 
stretching modes are Ub = 576.7 cm~^, Ug = 1015.5 cm~^, and coa = 1118.7 cm~^, respec- 
tively. The vibrational ground state (0,0,0) lies at 64220 cm~^ but remains invisible in a 
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perpendicular transition in which only states with odd number of antisymmetric stretch 
quanta Va are excited. Most intense lines in the spectrum belong to one of the three main 
progressions, assigned by visual inspection of the corresponding wave functions, although 
excitations are easily predictable from the FC geometries. The state ^S^ is linear, the an- 
gular dependence of the TDM ^^ is weak, and the bending vibration remains unexcited. 
In contrast, the equilibrium bond lengths, Ri = R2 = 2.4 ao, are substantially extended 
as compared to the equilibrium in X at Ri = R2 = 2.20 Cq, with the implication that the 
spectrum is built on stretching excitations. The nodal patterns of the vibrational eigenfunc- 
tions clearly demonstrate that the three main progressions, {vg, 0, 1), {vg, 0, 3), and {vg, 0, 5), 
involve consecutive excitations of the symmetric stretch mode accompanied by one, three, or 
five quanta of the antisymmetric stretch. Maximum intensity in each progression is reached 
for the same Vg = 7, and the intensity distributions around the maximum are approximately 
bell shaped. 

IV. SUMMARY AND CONCLUDING REMARKS 

This paper analyzes the absorption spectrum of carbon dioxide in the wavelength range 
120 nm — 160 nm by means of quantum mechanical calculations using PESs of five singlet 
valence electronic states and their coordinate dependent TDM vectors with the ground 
electronic state, constructed as described in paper I. The main results are as follows: 

1. The UV spectrum of five valence states, calculated for T = 190 K via Boltzmann 
averaging of optical transitions from many initial rotational states, reproduces the 
experimental two-band spectral envelope, the positions of the absorption maxima, 
their FWHMs, peak intensities, and frequencies of diffuse structures in each band. 

2. The two absorption bands have been assigned at several levels of graining. Electronic 
assignment, reflecting the strong HT effect and influenced by the vibronic interactions 
at CIs, is given in the adiabatic representation (low energy band: bent states 2M'/1 W; 
high energy band: linear states 3^A'/3^A"; gap between the two bands: state 2^A"). 
The parallel and perpendicular components of the full spectrum are isolated. The 
parallel transition, the RT coupling in which is suppressed for excitations from the 
ground vibrational state in X, is shown to be dominant for all but the highest photon 
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energies. Finally, individual diffuse structures are assigned using multicomponent wave 
functions of metastable resonance states. 

3. In the high energy band, the main progressions correspond to consecutive excitations 
of the pseudorotational motion around the closed loop of the CI seam; progressions 
differ in the number of nodes along the radial mode perpendicular to the closed seam. 
Systematic frequency growth in the pseudorotational progressions is interpreted as a 
spectroscopic hallmark of the closed seam. 

4. In the low energy band, the diffuse peaks, associated with bending excitations in the 
states 2^A' and l^A", are analyzed in both parallel and perpendicular transitions. 
Calculations reveal a dramatic difference in the spectra of the two states. The per- 
pendicular band (mainly iM"), bending progressions in which are nearly equidistant, 
is rather regular. In contrast, bending excitations in the parallel band (mainly 2^A'), 
are very anharmonic and irregular. This difference is interpreted as a manifestation 
of the carbene-type 'cyclic' OCO minimum in the PES of the 2^A' state. 

The present study focuses on spectroscopic aspects of the UV photodynamics of CO2. 
Mechanisms of photodissociation are discussed only insofar as they become apparent in 
the analysis of the absorption spectrum. This piece will be completed in the discussion 
of photofragment distributions which will be published separately. However, preliminary 
conclusions can be deduced from the density distributions of resonance wave functions ^n 
calculated and discussed above. Of interest here are the 2A' components of ^„ as depicted 
for example in Fig. [9j Although useless in establishing the vibrational assignments, they 
rather clearly illustrate the dissociation path taken by the molecule once it arrives in a bent 
adiabatic state. Because of CIs in linear CO2, the inner regions of states 2^A' and iM" are 
separated from the dissociation asymptotes by barriers. On the other hand, the 2^A' and 
iM" potentials along aoco in the direction of global bent equilibria are barrierless. This 
directional anisotropy is reflected in the probability density in the right panel of Fig. [9] The 
wave function traces out a path leading from the FC region to smaller angles; hardly any 
probability is found along the linear dissociation path. Substantial extension of the CO 
bond, seen in the left panel of Fig. |9} occurs predominantly in bent molecular configurations 
with aoco < 140°. This shape is representative of the 2A' or lA" components of ^„ both 
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below the CI energy (where the bent states are directly accessed) and above the CI energy 
(where they become populated indirectly via the CI funnels at linearity). 
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FIG. 1: (a) Experimental absorption cross section of CO2 {T = 195 K; in cm^) as a function 
of photon energy (in lO^cm"^). Numbers on the cross section axis denote powers of ten. The 
photon energy range studied in this work is indicated with two dashed hues, (b) Absorption 
spectrum, calculated at T = 200 K (blue line, vertically shifted by 0.9 • 10~^^ cm^), compared with 
the experimental spectrum (black line). Electronic assignment in terms of adiabatic states is given. 
Experimental data are from Ref. 1171 Theoretical spectrum is shifted by 400 cm~^ to lower energies. 

FIG. 2: Cuts through the PESs of the adiabatic states along (a) OCO angle (fixed coordinates are 
Ri = 2.2 oo and R2 = 2.3 oq) and (b) one CO bond distance (fixed coordinates are -Rco = 2.3 oq 
and aoco = 175°). The thick black line in both panels is the ground electronic state X; excited 
A' and A" states are shown with thin black and brown lines, respectively. Red lines and arrows 
illustrate the non-vertical excitation process and the effect of the TDM on the initial vibrational 
wave function tpQ in X. ipQ is shown in both panels with a dashed line, and the shapes of the initial 
state in the optical transitions via y and x components of the TDM, /iyVo (a) and fixi^o (b) are 
shown with solid lines. Arrows, indicating excitations, pass through the maximum of tpo and the 
maxima of iJ.yijJo and ^xi^o- 

FIG. 3: Spatial components fix (a,d), /Xy (b,e), and /i^ (c,f) of the ab initio TDM vectors for the 
adiabatic states 2^A' (purple), 3^yl' (brown), l^A" (black), 2^A" (purple), and ?^A" (brown) from 
the ground state X^A' . A' states are shown with solid. A" states with open symbols. Dependences 
on the CO bond distance are shown in (a) — (c) [fixed coordinates are R2 = 2.2 ao and aoco = 179°] 
and dependences on the OCO angle are shown in (d) — (f) [fixed coordinates are Ri = 2.2 oq and 
i?2 = 2.3 ao]. TDMs are measured in atomic units. The coordinate axes are sketched in (b): y 
is directed along the molecular figure axis, and z is normal to y in the molecular plane, and x is 
normal to the molecular plane. 
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FIG. 4: (a) The absorption spectrum a{Epii\Ji = 0) for a single initial state Jj = (blue) compared 
with the spectrum of Ref. [T7| measured at T = 195 K (black). Thin sticks indicate the major 
experimental peaks separated into a strong (brown) and a weak (blue) progression in the high 
energy band, and a strong progression in the low energy band (black). Their positions are taken 
from Table X of Ref. [T5l augmented with data from Ref. |T71 Panels (b-d) show decomposition of the 
calculated spectrum and vibronic assignments: (b) o"(£'ph| Ji = 0) spectrum (thin gray line) broken 
down into two low resolution components stemming from excitations of the bent {2^A' , l^A" , thick 
black line) and linear {3^A' , 3^A" , thick blue line) adiabatic states. Two spectra at the bottom 
are the parallel (crn, thick green line) and the perpendicular {a±, thick red line) components of 
a{EpYi\Ji = 0); (c) Decomposition of cjii (green) in terms of resonance states of the Hamiltonian Eq. 
( |11[ ). Gray sticks mark resonance positions (every 3d calculated state is shown). Black line is the 
sum of Lorentzians in Eq. (16); (d) Vibronic assignment of an (green) in terms of the most intense 



resonances shown at the bottom of the panel with thin sticks. Vertical shift of the calculated 
spectra is 0.9 • 10^^^ cm^ in panels (a) and (b), and 0.6 • 10^^*^ cm^ in panels (c) and (d). 

FIG. 5: Left and right columns depict 3A' adiabatic components of the 3D wave functions of 
resonance states in progressions (^^,0,0) and (w0,0, 1), respectively. Shown in molecular coordi- 
nates is one particular contour |^(-Ri, R2, aoco)| = const. The contour is viewed along the aoco 
axis and appears 'projected' onto the (i?i,i?2) plane. Shading emphasizes the 3D character of the 
plots. Middle column depicts contour maps of the eigenstates V'2d(-Rii -^2) calculated in the 2D 
single state model. All states belong to progression (t',/,,0)2D; red contours represent positive, blue 
contours negative values of V'2D- In all panels, a 2D contour map of the 3 A' adiabatic potential for 
c^oco = 179° (gray lines) and the CI seam (black thick line) are sketched. Axis tic labels are in 
aQ. States marked (u^, 0)20 are mixed with states {v^ — 4, 1)2^) shown in Fig. ^ The top middle 
panel shows author's drawing of a Mongolian hat. 

FIG. 6: Upper panels: 3 A' adiabatic components of the 3D wave functions of resonance states in 
the progression (^^,1,0). The state (3,1,0) could not be found. Lower panels depict contour 
maps of the eigenstates V'2d(-Ri;-R2) calculated in the 2D single state model and belonging to the 
progression (f^, 1)2D- The layout of all wave functions is the same as in Fig. I5| The 2D states are 
mixed with states {v^ + 4, 0)20; counterparts of those marked {v^, 1)20 are shown in Fig. [SJ 
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FIG. 7: Upper panels: Frequencies (a) in the major exprimental progression, (b) in the progression 
(f(^,0, 0) in the full calculations, and (c) in the progression (w0,O)2D in the 2D model. In (b), 
frequencies are lifted by 50 cm^^; in (c), the frequency scale is omitted. Arrows emphasize the 
positions of dips in the progressions. Lower panels: Potential of the 3^^' state (d) at aoco = 179° 
along the ab initio CI seam and (e) across the Mongolian hat top. Angular coordinate (j) along 
the seam in (d) is chosen such that (j) = 180° for Ri = R2 = 2.24 oq; the 'radial' coordinate in (e) 
runs along the antisymmetric stretch i?_, and the origin i?_ = is the point Ri = R2 = 2.5 oq. 
Potential curves are rised through the zero point energies of 'missing' coordinates, huj^j^ and fujo^/2 
in (d) and (e), respectively. Vibrational ladder in the 2D model is indicated in (d) and (e). Energy 
spacings forming two dips in the progression in panel (c) are highlighted with gray and yellow. 

FIG. 8: Shown left of the vertical dashed line are frequencies in the pure progressions of eigenstates 
^n in which the 2^A' adiabatic component (purple, blue, and brown lines and symbols) or the l^A" 
adiabatic component (red lines and symbols) are most populated. Open black symbols are the 
two experimental vibrational progressions (triangles and squares), as well as the frequency shift 
between their band maxima (circles) taken from Table 3 of Ref. [121 Shown right of the vertical 
dashed line are energy intervals between adjacent diffuse peaks in the experimental spectrum of 
Ref. [T7| (black open circles) and between most intense resonance states (green solid circles) . Note 
that the frequency scale changes in the upper part of the figure. 

FIG. 9: 2A' adiabatic component of the resonance state with Eq = 8.3795 eV and Fq = 62 cm^^. 
Shown is one particular contour |^(i?i, i?2i ooco)| = const., with the same constant in both 
panels. The contour is viewed along aoco (left panel) and along the CO bond (right panel). 
Shading emphasizes the 3D character of the plots. The 2D contour maps of the adiabatic state 
2^A! in the (i?i,i?2) plane and in the (00005-^2) plane are sketched in the left and right panels, 
respectively. Thick black lines trace out the probability density buildup away from the linear FC 
region. 

FIG. 10: 2 A' adiabatic components of the 3D wave functions of states in the progression (u^, 0, 0) 
(upper panels) and (0, 0, Va) (lower panels). The layout of all wave functions is the same as in the 
left panel of Fig. [9| 
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FIG. 11: 2A' adiabatic components of the 3D wave functions of states in the progression (0, Wf,, 0). 
In the left lower frame, the ground vibrational state in the carbene-type OCO minimum is shown 
and labeled (0, 0, 0)c- The layout of all wave functions is the same as in the right panel of Fig. M 



FIG. 12: Absorption spectrum of the state ^S„ and assignments of the three progressions {vs, 0, Va) 
in which the most intense absorption lines are found. 
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